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ABSTRACT 

We investigate the use of spiral arm pitch angles as a probe of disk galaxy mass profiles. We confirm 
our previous result that spiral arm pitch angles (P) are well-correlated with the rate of shear (S) in 
disk galaxy rotation curves, by using a much larger sample (51 galaxies) than used previously (17 
galaxies). We use this correlation to argue that imaging data alone can provide a powerful probe of 
galactic mass distributions out to large lookback times. In contrast to previous work, we show that 
observed spiral arm pitch angles are similar when measured in the optical (at 0.4 /xm) and the near- 
infrared (at 2.1 /Ltm) with a mean difference of 2?3 ± 2?7. This is then used to strengthen the known 
correlation between P and S using B band images. We then use two example galaxies to demonstrate 
how an inferred shear rate coupled with a bulge-disk decomposition model and a TuUy-Fisher derived 
velocity normalization can be used to place constraints on a galaxy's baryon fraction and dark matter 
halo profile. We show that ESO 582-G12, a galaxy with a high shear rate (slightly declining rotation 
curve) at ^ 10 kpc, favors an adiabatically contracted halo, with high initial NEW concentration 
(cvir > 16) and a high fraction of halo baryons in the form of stars 15 — 40%). In contrast, IC 
2522 has a low shear rate (rising rotation curve) at ~ 10 kpc and favors non-adiabatically contracted 
models with low NEW concentrations (cvir ~ 2 — 8) and a low stellar baryon fraction < 10%. 
Subject headings: dark matter — galaxies: fundamental parameters — galaxies: halos — galaxies: 
kinematics and dynamics — galaxies: spiral — galaxies: structure 



1. INTRODUCTION 

The correlation found between spiral arm pitch angle 
and rotation curve shear rate (Seigar, Block & Puerari 
2004; Seigar et al. 2005) suggests that there is a link 
between the tightness of spiral structure and the central 
mass concentration in spiral galaxies. The shear rate, S, 
a dimensionless quantity, can be measured directly from 
rotation curves and is defined as follows, 
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where A is the first Oort constant, to is the angular ve- 
locity and V is the velocity at a radius R. The shear 
rate depends upon the shape of the rotation curve. For 
a rotation curve that remains flat S = 0.5, for a falling 
rotation curve S > 0.5 and for a continually rising ro- 
tation curve S < 0.5. As the shape of a rotation curve 
depends on the mass distribution, the shear rate at any 
given position depends upon the mass within that radius, 
or the central mass concentration. As a result, the spiral 
arm pitch angle is dependent upon the central mass con- 
centration, and this is consistent with the expectations 
of most spiral density wave models (e.g. Bertin et al. 
1989a, b; Bertin 1991, 1993, 1996; Bertin & Lin 1996; 
Fuchs 1991, 2000), although density wave models pre- 
dict that pitch angles also depend on stability (i.e. the 
Toomre Q-parameter). 

The modal theory of spiral structure (Bertin et al. 
1989a, b; Bertin 1991, 1993, 1996; Bertin & Lin 1996) 
predicts that the tightness of the arms comes from the 
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central mass concentration. Galaxies with higher cen- 
tral mass concentrations, i.e. higher overall densities (in- 
cluding dark matter) , are predicted to have more tightly 
wound spiral arms. The models of Fuchs (1991, 2000) re- 
sult in disks with rigidly rotating spiral modes, wherein 
bulges act as inner reflectors of waves or modes induced 
by the swing amplification method, thus leading to modal 
spiral waves, which form as a result of Toomre swing am- 
plification (Toomre 1981). Fuchs (2000) adopts a stellar- 
dynamical analogue (Julian & Toomre 1966) of the Gol- 
dreich & Lynden-Bell (1965) sheet, which describes the 
local dynamics of an unbounded patch of thin, differen- 
tially rotating stellar disk. He then increases Q in the 
Goldreich & Lynden-Bell (1965) sheet, which has the af- 
fect that the bulge acts as an inner boundary. The result 
is that, instead of shearing density waves (as in the un- 
bounded sheet), spiral modes appear. These models best 
show how central mass concentration correlates with spi- 
ral arm pitch angle. If the disk is very light (low a where 
a is the disk density) the mode can be very tight, and 
we are in the domain of small epicycles. Formally, if the 
stability parameter Q — cK/irGa is close to unity, the 
value of c must also be small, where c is the radial ve- 
locity dispersion and k is the epicyclic frequency. If we 
increase the relative mass of the disk, we find a trend 
towards more open structures. 

Recent observational studies of spiral structure have 
highlighted a difference in the morphologies seen in the 
optical compared with the near-infrared. The classifica- 
tion of galaxies by Hubble type (Hubble 1926) is per- 
formed in the optical regime, where dust extinction still 
has a large effect, and where the light is dominated by 
the young Population I stars. In the near- infrared dust 
extinction is minimized and the light is dominated by old 
and intermediate age stars (Rhoads 1998; Worthey 1994). 
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Previous work has found that optically determined Hub- 
ble type is not correlated with the near-infrared morphol- 
ogy. There is not a good correlation between Hubble type 
and near-infrared bulge-to-disk ratio or near-infrared spi- 
ral arm pitch angle (de Jong 1996; Seigar & James 1998a, 
b). However, even in the optical regime, the correla- 
tion between quantitative pitch angle and Hubble type 
is weak (Kennicutt 1981), even though the tightness of 
the spiral arm pattern, as judged by eye, is one of the 
defining Hubble morphological criteria. Furthermore, it 
has been shown that near-infrared morphologies of some 
spiral galaxies can be significantly different from their 
optical morphologies (Block & Wainscoat 1991; Block et 
al. 1994; Thornley 1996; Seigar & James 1998a, b; Seigar, 
Chorney & James 2003). Galaxies with flocculent spiral 
structure in the optical may present grand-design spiral 
structure in the near-infrared (Thornley 1996; Seigar et 
al. 2003). These results suggest that the optical mor- 
phology gives incomplete information about the nature 
of spiral structure in disk galaxies. This is the basis of 
the near-infared studies of Block et al. (1999) and Seigar 
et al. (2005). 

However, work by Eskridge et al. (2002) showed that 
when assigning a Hubble type to galaxies in several dif- 
ferent wavebands, from the optical to the near-infrared, a 
correlation exists between the optically determined Hub- 
ble type and that determined in the near-infrared. This 
suggests that, while there may be large small-scale differ- 
ences between optical and near-infrared morphologies of 
spiral structure, the large-scale differences are minimal. 
As a result, it may be possible that there is no signifi- 
cant difference between the large-scale spiral arm pitch 
angle, regardless of whether it is determined in the opti- 
cal or near-infrared for the regular Grand-Design spirals 
investigated here. It should be noted that for flocculent 
spiral structure it is difficult to determine a pitch an- 
gle for the two-armed component in the optical, while 
a TO = 2 Grand-Design structure may be seen in the 
near-infrared (e.g. Thornley 1996; Grosbol & Patsis 1998; 
Seigar, Chorney & James 2003). As a result it may be 
possible to use optical imaging to expand the sample of 
galaxies having quantitative measurements of spiral arm 
pitch angle and therefore strengthen the correlation be- 
tween spiral arm pitch angle and rotation curve shear 
rate reported by Seigar et al. (2005). 

In addition to providing an important constraint on 
models of spiral arm formation, the correlation between 
spiral arm pitch angle and shear rate opens up a fun- 
damentally new approach for probing mass distributu- 
ions in spiral galaxies. This approach relies on imaging 
data alone without the need for full rotation curve infor- 
mation. Specifically, the TuUy-Fisher relation (Tully & 
Fisher 1977) for spiral galaxies coupled with the shear 
rate - pitch angle relation can be used to determine a 
rotation curve normalization and slope. In this paper we 
explicitly demonstrate how, given a bulge-disk decompo- 
sition, a pitch angle determination, and a Tully-Fisher 
normalizaion, one can constrain galaxy mass distribu- 
tions, dark matter halo concentrations, and other galaxy 
formation parameters. We obtain our constraints within 
the context of the the standard framework of disk forma- 
tion put forth by Fall & Efstathiou (1980) and Blumen- 
thal et al. (1986). 

In principle, the technique we demonstrate here can 



be applied generally to a large sample of galaxies and to 
galaxies at high redshift, when spiral arms are detected. 
This is the first in a series of papers in which we use spi- 
ral arm pitch angles to determine the mass distribution 
in spiral galaxies. In this paper we determine mass con- 
centrations in a few disk galaxies from their spiral arm 
pitch angles (or shear rates) and their disk masses and 
scalelengths (determined via a bulge-disk decomposition 
technique). In future papers we will apply the outline 
methods in this paper to a large sample of nearby disk 
galaxies and to a large sample of galaxies at higher red- 
shift in order to determine the evolution in the central 
mass concentrations of disk galaxies as a function of look- 
back time. 

This paper is arranged as follows. In section 2 we de- 
scribe the data used and the methods for determining 
spiral arm pitch angles and rotation curve shear rates. 
In section 3 we compare pitch angles in the near-infrared 
and optical. In section 4 we discuss the relationship be- 
tween spiral arm pitch angle and rotation curve shear 
rate found by Seigar et al. (2005) and we use more data to 
strengthen this correlation. In section 5 we use a bulge- 
disk decomposition routine and an adiabatic infall model 
for disk galaxy formation to investigate the link between 
shear rate and central mass concentration and we sum- 
marize our findings in section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

We have made use of H and B band images of 57 
galaxies from the Ohio State University Bright Spiral 
Galaxy Survey (OSUBSGS; e.g. Eskridge et al. 2002). 
These galaxies were chosen to be as face-on as possible 
(with a ratio of the minor- to major-axis b/a > 0.5) so 
that a comparison of spiral arm morphology could be 
made in the optical and near-infrared. 

We also include B band images of 31 galaxies that also 
have measured rotation curves (Mathewson et al. 1992; 
Persic & Salucci 1995). We observed these galaxies as 
part of an ongoing survey of the 600 brightest south- 
ern hemisphere galaxies. The observations make use of 
the CCD direct imaging camera at the 2.5-m du Pont 
telescope at the Las Campanas Observatory in Chile. 
Integration times used were 2 x 360 seconds. For this 
work, we selected from the full imaging dataset a sam- 
ple of galaxies having a ratio of the minor- to major-axis 
b/a> 1/3, so that spiral arm morphology can be studied, 
and rotation curves are also available. Of these galaxies, 
6 also have Ks band data. The Kg band images were 
observed using the Wide-field Infrared Camera (WIRC), 
also at the 2.5-m du Pont telescope. A 5-point dither 
pattern was used with integration times of 3 x 60 sec- 
onds in each position. The total exposure was therefore 
900 seconds. 

We also include B and Ks band images of a further 
3 galaxies from our survey, which do not have rotation 
curves. The selection for these galaxies satisfies the same 
selection criteria that was used to select the OSUBSGS 
galaxies above. 

2.1. Measurement of spiral arm pitch angles 

Spiral arm pitch angles are measured using the 
same technique employed by Seigar et al. (2005). A 
two-dimensional fast-Fourier decomposition technique is 
used, which employs a program described by Schroder 
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TABLE 1 

B BAND AND H BAND SPIRAL ARM PITCH ANGLES FOR 57 GALAXIES FROM THE OSUBSGS. 
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Note. — Column 1 lists the galaxy names; Column 2 lists the optical Hubble type from de Vaucouleurs et al. (1991; hereafter RC3); 
Column 3 lists the ratio of the minor- to major-axis, b/a from NED; Column 4 lists the position angle from RC3; Column 5 lists the three 
radial ranges used for the FFT analysis; Column 6 lists the B band pitch angle, Pb and column 7 lists the H band pitch angle, Ph- 



et al. (1994). Logarithmic spirals are assumed in the The amplitude of each Fourier component is given by 

decomposition. The resulting pitch angles are listed in 

Tables 1, 2 and 3. , Ei=i ELi (Inr, 6*) cxp -[i(m^p Inr)] 

A{m,p) = — f -J , 

(2) 
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TABLE 2 

B BAND AND Ks BAND SPIRAL ARM PITCH ANGLES FOR 9 GALAXIES FROM OUR SURVEY. 
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32 


21-175; 16-180; 26-170 


13.8±0.2 


13.0±0.7 


NGC 2082 


SABc 


0.94 


60 


50-130; 45-135; 55-125 


27.6±0.5 


20.9±0.6 


NGC 2280 


SAcd 


0.49 


163 


55-117; 50-122; 60-112 


24.2±1.7 


22.4±2.0 


NGC 2417 


SABbc 


0.68 


81 


25-71; 20-76; 30-66 


24.0±0.7 


22.2±3.7 


NGC 2935 


SABb 


0.78 





43-93; 38-98; 48-88 


14.1±0.2 


16.5±1.1 


NGC 3318 


SABb 


0.54 


78 


23-89; 18-94; 28-84 


36.9±6.5 


36.9±2.2 


NGC 3450 


SBb 


0.88 


140 


28-48; 25-51; 31-45 


9.1±0.4 


14.5±0.3 


NGC 4050 


SBab 


0.68 


85 


70-91; 67-94; 73-88 


8.9±0.7 


8.8±1.6 



Note. — Column 1 lists the galaxy names; Column 2 lists the optical Hubble type from RC3; Column 3 lists the ratio of the minor- to 
major-axis, b/a from NED; Column 4 lists the position angle from RC3; Column 5 lists the three radial ranges used for the FFT analysis; 
Column 6 lists the B band pitch angle, Pg and column 7 lists the Ks band pitch angle, Px- 



where r and 9 are polar coordinates, I{hir,9) is the in- 
tensity at position (In 9) , m represents the number of 
arms or modes, and p is the variable associated vi^ith the 
pitch angle P, defined by tanP = —(m/p). Throughout 
this work we measure the pitch angle P of the m = 2 
component. The resulting pitch angle measured using 
equation El is in radians, and this is later converted to 
degrees for ease of perception. 

The range of radii over which the Fourier fits were ap- 
plied were selected to exclude the bulge or bar (where 
there is no information about the arms) and to extend 
out to the outer limits of the arms in our images, in 
such a way that the 10 kpc radius fell approximately 
in the middle of this range. The radial extent of the 
bar was measured manually (see e.g. Grosbol, Patsis & 
Pompei 2004), and the inner radial limit applied to the 
FFT was chosen to be outside this radius. Physical dis- 
tances are calculated using a Hubble constant Hq — 75 
km s~^ Mpc~^ and recessional velocities from the NASA 
Extragalactic Database (NED). Pitch angles are then de- 
termined from peaks in the Fourier spectra, as this is the 
most powerful method to find periodicity in a distribu- 
tion (Considere & Athanassoula 1988; Garcia-Gomez & 
Athanassoula 1993). The radial range over which the 
Fourier analysis was performed was chosen by eye and is 
probably the dominant source of error in the calculation 
of pitch angles, as spiral arms are only approximately 
logarithmic and sometimes abrupst changes can be seen 
in spiral arm pitch angles (e.g. Seigar & James 1998b) As 
a result, three radial ranges were chosen for each galaxy, 
and a mean pitch angle and standard error calculated for 
every object. 

The images were first projected to face-on. Mean un- 
certainties of position angle and inclination as a function 
of inclination were discussed by Considere & Athanas- 
soula (1988). For a galaxy with low inclination, there 
are clearly greater uncertainties in assigning both a posi- 
tion angle and an accurate inclination. These uncertain- 
ties are discussed by Block et al. (1999) and Seigar et al. 
(2005), who take a galaxy with low inclination (< 30°) 
and one with high inclination (> 60°) and varied the in- 
clination angle used in the correction to face-on. They 
found that for the galaxy with low inclination, the mea- 



sured pitch angle remained the same. However, the mea- 
sured pitch angle for the galaxy with high inclination var- 
ied by ±10%. Since inclination corrections are likely to 
be largest for galaxies with the highest inclinations cases 
where inclination is > 60° are taken as the worst case 
scenario. For galaxies with inclination i > 60° we take 
into account this uncertainty. Our deprojection method 
assumes that spiral galaxy disks are intrinsically circular 
and flat in nature. 

2.2. Measurement of shear rates 

31 of the galaxies observed here have Ha rotation curve 
data measured by Mathewson ct al. (1992) and Persic & 
Salucci (1995). These rotation curves are of good quality 
with an rms error < 10 km s~^, and an error associated 
with folding the two sides of the galaxy also < 10 km 
s~^. These rotation curves have been used to estimate 
the shear rates in these galaxies, using the same method 
used by other authors (e.g. Block et al. 1999; Seigar et 
al. 2005; Seigar 2005). 

Using equation n we have calculated the shear rates 
for these galaxies, over the same radial ranges for which 
the Fourier analysis was performed and pitch angles cal- 
culated. We have selected several different radial ranges, 
just as in the Fourier analysis, and we present mean shear 
rates and standard errors. The dominant sources of er- 
ror on the shear rate are the rms error in the rotation 
curve and the error associated with folding the two sides 
of the galaxy. This is typically < 10%. In order to calcu- 
late the shear rate, the mean value of dV/dR measured 
in km arcsec"^ is calculated by fitting a line of con- 
stant gradient to the outer part of the rotation curve 
(i.e. past the radius of turnover and any bar or bulge 
that may exist in the galaxy). Mean shear rates are 
then calculated from shear rates measured over three ra- 
dial ranges, corresponding to the same radial ranges over 
which the Fourier analysis was performed. The resulting 
shear rates are listed in Table 3. 

The choice of using 10 kpc as the radius at which to 
measure the shear and spiral arm pitch angle is somewhat 
arbitrary. The correlation between these two quantities 
(Seigar et al. 2005 and Figure 3) only exists when they 
are measured at a physical radius, chosen indepently of 
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B band pitch angle in degrees 

Fig. 1. — A 1:1 correlation between optical spiral arm pitch angle 
as measured in the B band and near-infrared spiral arm pitch an- 
gle as measured in the H band for OSUBSGS galaxies and the Ks 
band for our galaxies. This correlation shows that pitch angles are 
similar whether measured in the near-infrared or optical regime. 
Only one object shows a difference of > 10°. The solid squares 
represent galaxies from the OSUBSGS and the open squares rep- 
resent galaxies from our survey. The solid line represents where a 
1:1 correlation lies. The dashed line is the best-fit line to the data. 



the disk scalelength. If we measure pitch angles and shear 
at a radius chosen using the disk scalelength, the correla- 
tion no longer exists, and the range of shear is narrower 
(between 0.4 and 0.6). As a result, we chose to measure 
shear and pitch angles at a radius independent of the disk 
scalelength. We chose 10 kpc (although it could just as 
easily be 8 kpc or 12 kpc). 

The fact that this correlation only exists for a physical 
radius, suggests that other quantities may be important 
in the determination of shear and pitch angle (e.g. the 
disk scalelength itself or the stability of the disk as a 
function of radius). Indeed, when modeling galaxies we 
use a 1-dimensional bulge-disk decomposition and use 
the derived disk scalelength in our codes. 

3. A COMPARISON OF OPTICAL AND NEAR-INFRARED 
PITCH ANGLES 

Figure ^ shows of plot of B band pitch angle versus 
near-infrared pitch angle for the 57 face-on and nearly 
face-on galaxies taken from the OSUBSGS and 9 face-on 
and nearly face-on galaxies taken from our survey. This 
plot shows a very tight 1:1 correlation between the optical 
and near-infrared spiral arm pitch angles (correlation co- 
eSicient=0.95; significance>99.99%), showing that pitch 
angles are similar whether measured in the optical or 
near- infrared regimes. In this plot, only one galaxy has 
a pitch angle that shows a difference of > 10° when mea- 
sured at B compared with its near-infrared measurement. 
The outlying galaxy is NGC 613, which appears to have 
a strong 4-armed component in the optical, whereas this 
is a much weaker feature in the near-infrared. 

This result is in contrast to the claims of Block et 
al. (1999) and Block & Puerari (1999), who argue that 
there is usually a large difference between the spiral arm 
pitch angle when measured in the optical and the near- 
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Difference between optical and near-IR pitch angles 

Fig. 2. — Histogram showing the distribution in the absolute 
difference measured between the B band and H band pitch angles. 



infrared. Our result is consistent with the result that 
Hubble types assigned to galaxies in the optical and the 
near-infrared correlate well with each other (Eskridge et 
al. 2002). 

We now analyze the absolute difference between the 
pitch angles measured the the B band and those mea- 
sured in the near-infrared. The mean of this difference 
in pitch angles is 2?3 ± 2?7, where the error is a standard 
deviation, and this is therefore consistent with the same 
pitch angle being measured, regardless of the waveband. 
The mode of the distribution is in the range 0.5—1.0 and 
the median is 1.46. Furthermore a linear fit to the data 
in Fig n yields the following relationship, 

Pnir = (0.50 ± 1.25) + (0.98 ± 0.04)Popt, (3) 

where Pnir is the pitch angle measured in the near- 
infrared and Popt is the pitch angle measured in the op- 
tical. The point at which this line of best fit intersects 
the P/v/i?-axis is consistent with zero, and the gradient 
of the line is consistent with unity, and so the line of 
best fit is consistent with a 1:1 correlation between spi- 
ral arm pitch angles in the optical and near-infrared. It 
should be noted that although the difference between 
optical and near-infrared spiral arm pitch angles is con- 
sistent with zero from a statistical analysis, minor sys- 
tematic differences cannot be excluded. Small differences 
between optical and near-infrared pitch angles have been 
measured of 1° - 2° (e.g. Grosbol & Patsis 1998). Fur- 
thermore, color gradients across arms may be caused by 
such small differences in pitch angles (e.g. Gonzalez & 
Graham 1996). 

Of course, it is still a well-known fact that spiral 
structure measured in the near-infrared appears much 
smoother than that measured in the optical (e.g. Block & 
Wainscoat 1991; Block et al. 1994; Thornley 1996; Seigar 
et al. 2003). Even in galaxies that appear grand-design in 
the optical, small spokes and bifurcations are clearly seen 
along the length of the dominant arms. These spokes 
and bifurcations all but disappear in the near-infrared. 
However, our result suggests that the pitch angle of the 
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Shear rate 



Fig. 3. — Spiral arm pitch angle versus rotation curve shear rate, 
showing a strong correlation. The solid squares represent galaxies 
with data measured by Block et al. (1999), the open squares are 
galaxies from Seigar et al. (2005), and the open triangles represent 
the data from the present sample. 



scheme introduced by Block & Puerari (1999). These 
authors argue that their dust-penetrated classification 
scheme is needed, because spiral arm pitch angles are 
different in the near-infrared when compared to the op- 
tical, and as a result they classify galaxies into three bins, 
based on their pitch angles when measured in the near- 
infrared. However, we have shown (in Figure ^ that 
spiral galaxies with a relatively strong two-armed spiral 
pattern, have optical and near-infrared pitch angles that 
are nearly always very similar. This, therefore, seems 
to lessen the need for the dust-penetrated classification 
scheme. In Figure 01 however, we show that a correla- 
tion still exists between shear rate and spiral arm pitch 
angle, no matter what waveband the pitch angle is mea- 
sured in. This suggests that classifying galaxies based on 
their pitch angles has some physical basis, because it tells 
us something about their central mass concentrations. 

Given the correlation between shear rate and pitch 
angle we derive the following expression relating these 
terms 

P = (64.25 ± 2.87) - (73.24 ± 5.53)5', (4) 

where P is the pitch angle in degrees and S is the shear 
rate. 



large-scale spiral structure remains constant with wave- 
length in the range 0.4/i,m< A <2.2/im for spirals with a 
relatively strong two-armed pattern in the optical. For 
weaker, more flocculent galaxies there may well be sig- 
nificant differences due to dust. This has serious impli- 
cations for the physical meaning of the dust-penetrated 
classification scheme (Block & Puerari 1999; Block et al. 
1999). These implications are discussed in section 0] 

4. THE RELATIONSHIP BETWEEN SHEAR RATE AND 
PITCH ANGLE 

The correlation between spiral arm pitch angle and ro- 
tation curve shear rate presented by Seigar et al. (2005) 
only had 17 data points, including 3 from an earlier study 
by Block et al. (1999). Here we include a further 31 
galaxies with measure pitch angles and shear rates. Fig- 
ure O shows a plot of shear rate versus spiral arm pitch 
angle. A very good correlation still exists (correlation 
cocfficient=0.89; significance=99.75%). 

The shape of a rotation curve, beyond the turnover 
radius, is determined by the amount and distribution of 
matter contained in a galaxy. The increase of shear rates 
from low to high dictates a change in mass distributions 
from small to large central mass concentrations. The cor- 
relation with pitch angle is therefore interpreted as fol- 
lows: galaxies with higher rates of shear present a larger 
central mass concentration and more tightly wound arms. 
In contrast, open arm morphologies are associated with 
lower rates of galactic shear and lower central mass con- 
centrations. This correlation was alluded to by the pio- 
neering work of Lin & Shu (1964) and by the later spec- 
troscopic study of Burstein & Rubin (1985). It is in 
agreement with modal theories of spiral structure (e.g. 
Bertin et al. 1989a, b; Berlin & Lin 1996) and other nu- 
merical models based on the modal theory (e.g. Fuchs 
1991, 2000). 

The correlations shown in Figures Q and |31 have signif- 
icant implications for the dust-penetrated classification 



5. CONNECTING PITCH ANGLES TO CENTRAL MASS 
CONCENTRATIONS 

Our goal now is to demonstrate how spiral arm pitch 
angle measurements coupled with a bulge-disk decompo- 
sition can be used to constrain the overall mass distri- 
bution in spiral galaxies. We outline the method for two 
example galaxies with pseudo-bulges (see Kormendy & 
Kennicutt 2004 for a review about pseudo-bulges). The 
two galaxies we chose are chosen such that one galaxy, 
ESO 582-G12, has a flat rotation curve, with a shear, 
S = 0.52 ±0.05 in the middle of the range, and the other 
galaxy, IC 2522, has a continually rising rotation curve, 
with a shear, S — 0.39 ±0.01, towards the low end of the 
range of shears. As a normalization for each galaxy we 
take its measured rotation speed at 2.2 disk scale lengths, 
V2.2 (listed in the last column of Table 4), and include a 
bulge-disk decomposition in order to estimate the bary- 
onic contribution to each rotation curve. In future work 
paper we plan to use the same methodology for a large 
sample of galaxies. 

For each of our galaxies we produce surface bright- 
ness profiles using the IRAF Ellipse routine, which 
fits ellipses to isophotes in an image, using an itera- 
tive method described by Jedrzejewski (1987). From 
the surface brightness profile we then determine the disk 
and bulge B-band luminosity using an exponential disk 
and an exponential (pseudo) bulge. In future papers a 
more general Sersic parameterization will be used for the 
bulge component. We utilize a 1-D bulge-disk decom- 
postion routine, which performs Levenberg-Marquardt 
least-squares minimization. Explicity we fit an exponen- 
tial profile for each bulge via 

fl{R) = fie cxp (-1.679[(i?/i?e) - 1]), (5) 

where Re is the effective radius, containing 50% of the 
total light of the bulge and fj,e is the surface brightness 
at Re- We fit the disk component using 

/i(i?) = ^0 exp {-R/h), (6) 
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TABLE 3 

Pitch angles and shear rates for 31 galaxies from our survey. 



Galaxy 


Hubble 


b/a 


PA 


Radial Ranges 


Pitch 


Shear 


Name 


Type 


(degrees) 


(arcsec) 


Angle 
(degrees) 


Rate 




SAc 


0.76 


171 


97—1 01 ■ 29 1 nfi- '^9— Qfi 


23.7ibl.l 




PSD 1 91 n9R 


SBc 


0.63 


115 


o± oOj 'lu yu, ou 




U. ( U^U. U1 




SAc 


0.64 


48 


32—1 02- 27-1 07- 37—97 


.U_1_U. u 


n ^o+n 05 

U. OiJ _1_U. UO 


TP 9^99 


SAcd 


0.71 


Q 


9'^— 77- 90—89- '^O— 79 


OO. ± . u 


n ^Q-i-n 01 

u. oy ^u. UJ- 


TC 2^37 


SABc 


0.65 


26 


"^0— 89- 9^-87- ^^^-77 

ijyj Oil) 1 ^ tjo 1 ( 


9.3di0.3 


n fiQ+0 09 

U. UO_l_U. Uii 




SAc 


0.39 


23 


99_Q9- 1 7_Q7- 97—87 




n '^fi-i-0 0"^ 

U. OU^U. UO 


TP 4.^'^^^ 


SABc 


0.77 


50 


A9_fiR- Q7_71 - 47_R1 


OO.O^O . o 


n ^f^-l-0 09 
u. oo^u. uz 


TP zL^^O)^ 


SAc 


0.42 


45 


1 1 (1— '^n- 9(1— zLfl 


14.1ib0.4 


n fic;-i-n 09 

u. uo^u. uz 


NPP 1 '^n 


SBbc 


0.50 


118 


70— 19fi- fic;_1 m ■ 7f;_191 


8.4di0.1 


n fic;-i-n O"^ 

U. UO^U. UO 


NPP 1 '^l 


SBbc 


0.42 


75 


'^0— ^9- 97— '\'^—AQ 


36.1ibl.5 


n '^9-1-0 09 


NTPP f^7R 


SABc 


u.oo 


1 1 n 


fif^ 19fi- fin 1*^0- 70 190 

IZu, UU lOU, lU J-ZU 


1 8 O-kO 9 
-LO.U^U. z 


n fi4-i-n ofi 

U-U4^U-UU 


NPP Qn^^ 


SAc 


0.43 


75 


•^fi- 1 i^O- '^1— fil— Id'^ 


12.9ib0.4 


u. oy^u. U4 


NPP 1 9'^9 


SABc 


0.88 


108 


^'i—^A^■ a8— izLfi- '^8— i'?fi 


19.3ib2.0 


n fifi-l-n Ofi 

U. UU^U. UU 


NCiC 1 2q2 


SAc 


0.43 




c;Q_i7'i- A8— 1 78- '^8— 1 fi8 

OO ± ( Oj 4:0 X 1 O, OO -1-UO 


29.8ibl 


n 48-1-0 04 

U.40^U. U4 


NPP 1 '^nn 


SBbc 


0.66 


106 


80— llfi- 7^1— 191- S'^— 1 1 1 

ou J-±U, i O ±Z±, OO ±-L± 


31.7ibl.l 


n f^O-l-n 0"^ 

U. OU^U. UO 


IVPP 1 '^'^'^ 


SAbc 


41 


-LOO 


f\'\ A'^Q- AAA- fi8 1 
UO lOy, OO ±44:, UO 104 


■^fi fi-H 
OO.U^ 1 .u 


n '^4-i-n 0^^ 

U.04^U.U0 


NPP 1 '^fi^i 


SBb 


0.55 


32 


fii_i9Q- f^fi—l '^4- fifi— 1 94 

U J- -L y ( ou ±04, UU J. ^4 


35.4ibl.7 


n f^'^-i-o 0*^ 


NOP 1 ^^Q 


SBcd 


0.57 


64 


82-156- 77-161- 87-151 

Oil -LiJU) 1 1 -LU-L, Ol -L O i 


20.4ib0.4 


f) fi8-i-n ofi 

U. UO_l_U. UU 


NPP 1 '^fifi 




0.80 


60 


7'i_lQ'i- fi8— 1 '^8- 78—198 
J-OOj UO ±00, (O ±zo 


n+O '\ 
ou.u^u. o 


n f^o-i-n 04 

U. OU^U. U4 


NGC 1964 


SABb 


0.38 


32 


21-175; 16-180; 26-170 


13.8±0.2 


0.61±0.02 


NGC 2280 


SAcd 


0.50 


163 


55-117; 50-122; 60-112 


24.2±1.7 


0.32±0.05 


NGC 2417 


SABbc 


0.68 


81 


25-71; 20-76; 30-66 


24.0±0.7 


0.63±0.03 


NGC 2835 


SABc 


0.67 


8 


110-240; 105-245; 115-235 


19.5±0.7 


0.62±0.06 


NGC 2935 


SABb 


0.78 





43-93; 38-98; 48-88 


14.1±0.2 


0.73±0.03 


NGC 3052 


SABc 


0.65 


102 


25-57; 20-62; 30-52 


19.8±0.7 


0.63±0.02 


NGC 3054 


SABbc 


0.61 


118 


44-84; 39-89; 49-79 


42.9±3.9 


0.37±0.02 


NGC 3223 


SAbc 


0.61 


135 


32-74; 27-79; 37-69 


10.7±2.0 


0.72±0.02 


NGC 3318 


SABb 


0.54 


78 


23-89; 18-94; 28-84 


36.9±6.5 


0.53±0.03 


NGC 5967 


SABc 


0.59 


90 


46-76; 41-81; 51-71 


47.3±0.5 


0.27±0.02 


NGC 7083 


SABc 


0.59 


5 


32-78; 29-81; 35-75 


26.7±1.4 


0.40±0.04 


NGC 7392 


SBab 


0.62 


123 


30-66: 25-71; 35-61 


24.6±2.0 


0.54±0.02 



Note. — Column 1 lists the galaxy names; Column 2 lists their Hubble typos from R,C3; Column 3 lists the ratio of the minor- to 
major-axis from NED; Column 4 lists the position angle from RC3; Column 5 lists the three radial ranges used for the EFT analysis; 
Column 6 lists the pitch angle and Column 7 lists the shear rate. 



where /io is the central surface brightness and h is the 
scalelength of the disk. The results of this bulge-disk 
decompostion can be seen in Figure^ and Table 4. 

We assign masses to the disk and bulge components 
using a range of stellar mass-to-light ratios from Bell 
& de Jong (2001). Specifically, in our rotation curve 
models we allow mass-to-light ratios of (M/L) = 1.0, 1.3 
and 1.6 (measured in B band solar units), and use 
our photometrically-derived disk and bulge light profiles 
Lb = Ldisk + Lhuige to determine the stellar mass con- 
tributuion to each rotation curve: — {M/L)Lb- 

We now explore a range of allowed dark matter halo 
masses and density profiles, adopting two extreme mod- 
els for disk galaxy formation. In the first, we assume 
that the dark matter halos surrounding these galaxies 
have not responded significantly to the formation of the 
disks, i.e. adiabatic contraction (AC) does not occur. We 
refer to this as our "non-AC" model. In this case, the 
dark matter contributuion to each galaxy rotation curve 
is described by a density profiles that mirrors those found 
in dissipationless dark matter simulations: 



p{r) 



Ps 



n2 ' 



(7) 



{r/rs) (1 + r/vs) 

where Tg is a characteristic "inner" radius, and ps a cor- 
responding inner density. Here we have adopted the pro- 



file shape of Navarro, Frenk & White (1996; hereafter 
NEW). The NEW profile is a two-parameter function 
and is completely specified by choosing two independent 
paramters, e.g. the virial mass Mvir (or virial radius 
i?vir) and concentration Cvir = Rvh/fs define the profile 
completely (see Bullock et al. 2001b, for a discussion). 
Similarly, given a virial mass Mvir and the dark matter 
circular velocity at any radius, the halo concentration 
Cvir is completely determined. 

In the second class of models we adopt the scenario of 
"adiabatic contraction" (AC) discussed by Blumcnthal 
et al. (1986; see also Bullock et al. 2001a, Pizagno et al. 
2005). Here we assume that the baryons and dark mat- 
ter initially follow an NEW profile and that the baryons 
cool and settle into the halo center slowly compared to 
a typical orbital time. This slow infall provokes an adi- 
abatic contraction in the halo density distributuion and 
gives rise to a more concentrated dark matter profile. 
The idea of adiabatic contraction was originally discussed 
as to explain the "conspiracy" between dark halos and 
disk sizes which gives rise to a featureless rotation curve 
(Rubin et al. 1985) but has since proven to be remark- 
ably accurate in describing the formation of disk galaxies 
in numerical simulations (e.g. Gnedin et al. 2004, and 
references therein), although the degree to which this 
process operates in the real universe is currently uncer- 
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Radius (arcsec) Radius (arcsec) 

Fig. 4. — Bulge-disk decomposition for two example spiral galaxies as observed in the B-band. The top row shows a B band image of 
each galaxy and the bottom row shows the surface brightness profile and results of the decomposition. The effective radii Re and disk 
scalelengths h arc measured in arcsec. The dashed line is the best fit bulge model, the dotted line is the best fit disk model and the solid 
line is the best combined bulge+disk model. Left panel: ESO 582-G12. Right Panel: IC 2522. The images are aligned such that North is 
up and the field is 4.5x4.5 arcmin. 



tain. For example, Button et al. (2005), showed that 
adiabatic contraction models are inconsistent with the 
rotation curves measured, and the expected NFW con- 
centrations, for a sample of 6 galaxies. They suggest that 
mechanisms such as stellar feedback and stellar bars may 
result in less concentrated halos than predicted by adia- 
batic concentration. 

In our AC model, we take the contraction into account 
following the prescription of Blumenthal et al. (1986). 
Note that Gnedin et al. (2004) advocate a slightly mod- 
ified prescription, but the differences between the two 
methods arc small compared to the differences between 
our "AC" model and our non-AC model. In principle, 
any observational probe that can distinguish between AC 
and non-AC type scenarios provides an important con- 
straint on the nature of gas infall into galaxies (i.e. was 
it fast or was it slow?). 



For each galaxy we iterate over the central and ±1 — cr 
values found in the bulge-disk decompositions for h and 
Ldisk and explore the three values of mass-to-light ratio 
discussed above [M/L) = 1.0, 1.3, and 1.6. In each case 
we assume that M/L remains constant with radius. For 
each choice of bulge-disk model parameters and mass-to- 
light ratios we allow the (initial) halo NFW concentration 
parameter to vary over the range of viable values: Cyir = 
3-31 (Bullock et al. 2001b). We then determine the 
halo virial mass M^-„ necessary to reproduce the value 
of V2.2 for the galaxy and determine the implied fraction 
of the mass in the system in the form of stars compared 
to that "expected" from the Universal baryon fraction: 
= M*/(/6Mvir). We make the (rather loose) demand 
that /* lies within the range of plausible values /* = 
0.01 - l/b. 

For each chosen value of Cyir and adopted disk forma- 
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TABLE 4 

Galaxy parameters used in model constraints 



Galaxy Name 


Distance 


h 


h 


B/D 




V2.2 




(Mpc) 


(arcsec) 


(kpc) 




(xIOIOLq) 


(km s^-"-) 


ESQ 582-G12 


31.0±0.1 


36.5±3.8 


5.48±0.57 


0.106 


1.27±0.11 


145 


IC 2522 


40.3±0.1 


20.4±2.1 


3.98±0.41 


0.164 


1.55±0.14 


216 



Note. — Column 1 lists the galaxy names; Column 2 lists the distances to the galaxies in Mpc, calculated using a Hubble constant 
Hq = 75 km s~^ Mpc~^; Column 3 lists the disk scalelengths, h, in arcsec; Column 4 lists the disk scalelengths, h, in kpc; Column 5 
lists the bulge-to-disk ratio , B/D; Column 6 lists the disk B-band luminosity, Ltjisfc; a^nd Column 7 lists the rotational velocity at 2.2 
disk scalelength, V2.2. 
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Fig. 5. — Model Results for ESO 582-G12. Left: NEW concentration versus shear, center, stellar baryon fraction versus shear and right: 
virial mass versus shear. The open squares represent the model without an adiabatically contracted dark matter halo (non-AC) and the 
solid circles represent the model with adiabatic contraction (AC). The vertical lines present the measured shear (solid line) with the l-cr 
error added and subtracted (dashed lines). 



tion scenario (AC or non-AC) the chosen V2.2 constraint 
defines the rotation curve completely, and thus provides 
an implied shear rate at every radius. The three panels of 
Figures 13 and El show the results of this exercise for ESO 
582-G12 and IC 2522 respectively. In each panel, open 
symbols are for the non-AC model and the solid sym- 
bols reflect the AC assumption. Each point represents a 
distinct input combination of h, Ldisk, and (M/L). The 
measured shear rate is illustrated by a solid vertical line 
and the ±1 — cr range in the observe shear rate for each 
galaxy is shown by the two vertical dashed lines in each 
panel. 

Consider first the left hand panel of Figure |H Here 
we plot the dark matter halo concentration parameter 
versus the shear rate measured at 10 kpc. More con- 
centrated halos generally produce higher shear rates, as 
expected. It can be seen that for a given NFW concentra- 
tion, Cyir, several values of shear arc possible. This is due 
to changes in the baryon contrbution (i.e. the disk mass 
and disk scalelength) to the rotation curve. Whether an 
increase in the baryon contribution causes the shear to 
increase or decrease depends on the size of the disk (i.e. 
the disk scalelength). The same is true for Figure 6. In 
the AC model (solid), Cvir refers to the halo concentra- 
tion before the halo is adiabatically contracted. This is 
why the AC points tend to have higher shear values at 
fixed Cvir compared to the non-AC case. Every point 
(or h, Ldisk, {M / L) input combination) has an associ- 
ated stellar baryon fraction and dark halo virial mass 
Mvir- These values are plotted versus shear rate in the 



middle and right-hand panels in Figure [S] Observe that 
with only the V2.2 constraint imposed (all points), then 
a wide range of dark matter halo properties are allowed. 
Once we constrain the models by forcing the predicted 
shear to be consistent the observed range S — 0.47 — 0.57, 
we favor models that include adiabatic contraction. The 
dark matter of halo of ESO 582-G12 is constrained to be 
of relatively low mass, Mvir — 5 — 8 x lO^-'^M©, and to 
have a fairly high NFW concentration, Cvir >~ 16. This 
is consistent with the high-end of the the Cvir distribution 
predicted for LCDM l|Bullock et al. 1120019 ). The stellar 
mass in the galaxy is between ~ 15 and 40% of the the 
baryonic mass associated with its dark matter halo. 

Figure El shows analogous results for IC 2522. Un- 
like ESO 582 G12, this galaxy prefers a more massive, 
low-concentration halo: Mvir — 1.5 — 6 x lO^^M© with 
Cvir — 2 — 8, at the low end of the Cvir distribution 
l|Bullock et al. ll2001bD . The imphed stellar mass is quite 
small: < 10% of the universal baryon budget has ended 
up on stars. Furthermore, the measured shear rate of 
this galaxy favors our non-AC model. This is in quali- 
tative agreement with results presented by Dutton et al. 
(2005). 

Model results for the central mass concentration in 
each galaxy are shown in Figures and |H1 Here Ctot 
and Cdm correspond to the fraction of the total and dark 
matter mass contained within central 10 kpc of each sys- 
tem. We see that ESO 582-G12 contains roughly 20% 
of its mass within 10 kpc, while IC 2552 is quite diffuse 
with Ctot < 5%. 
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Fig. 6. — Model Results for IC 2522. Left: NFW concentration versus shear, center: stellar baryon fraction versus shear and right: virial 
mass versus shear. The open squares represent the model without an adiabatically contracted dark matter halo (non-AC) and the solid 
circles represnt the model with adiabatic contraction (AC). The vertical lines present the measured shear (solid line) with the 1-cr error 
added and subtracted (dashed lines). 
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Fig. 7. — Model results for ESO 582-G12: Left: Central mass concentration (mass fraction within 10 kpc) versus shear, right: dark 
matter concentration (dark matter mass fraction within 10 kpc) versus shear. The open squares represent the model without adiabatic 
contraction of the stellar halo and solid circles represnt the model with adiabatic contraction. The vertical lines present the measured shear 
(solid line) with the l-cr error added and subtracted (dashed lines). 



What drives these differences? Take IC 2522 for ex- 
ample, which has a very low concentration. This arises 
because of the combination of its (small) disk scalelength 
and the fact that the rotation curve is rising at 10 kpc 
(implied by a shear of 0.39±0.1). The rotation curve of 
the stellar disk alone peaks at 2.2 scalelengths (~ 8.8 
kpc), and will be falling at 10 kpc (where the shear is 
measured). The measured shear of 0.39±0.01 at 10 kpc 
implies that the total rotation curve is rising at 10 kpc, 
which means the dark matter must be very extended (or 
of low concentration). ESO 582-G12, on the other hand, 
has a falling rotation curve at 10 kpc, but a larger disk 
scalelength. This requires a more concentrated halo to 
explain the observed shear rate. 

The two galaxies chosen for this modeling are, in many 
ways, quite comparable. Their Hubble types, B/D ra- 
tios and disk luminosities are very similar. The only 
significant difference is in their measured rates of shear. 



The differences in their favored formation models (AC for 
ESO 582-G12 versus non-AC for IC 2522), halo masses 
and concentrations, comes down to a different measure- 
ment for their shear and the sizes of their disks. It is 
quite interesting that two galaxies, which seem so simi- 
lar in size and Hubble type, seem to inhabit very differ- 
ent types of dark matter halos. It is tempting to spec- 
ulate that the shear rate itself provides an underlying 
physical driver to push galaxies in drastically different 
dark matter halos towards similar luminosities. How- 
ever, the problem would then be what stops such a drive 
and makes the galaxy luminosities so similar, when the 
halo masses show large differences. 

This clearly demonstrates that the shear rate adds an 
important constraint on galaxy formation models com- 
pared to what can be learned from standard Tully-Fisher 
constraints alone. Indeed, the results for IC 2522 and 
ESO 582-G12 are tantalizing. The first, with a low shear 
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Fic;. 8. — Model results for IC 2522. Left: Central mass concentration (mass fraction within 10 kpc) versus shear, right: dark matter 
concentration (dark matter mass fraction within 10 kpc) versus shear. The open squares represent the model without adiabatic contraction 
of the stellar halo and solid circles represnt the model with adiabatic contraction. The vertical lines present the measured shear (solid line) 
with the 1-cr error added and subtracted (dashed lines). 
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Fig. 9. — Observed rotation curves with the overlaid model rotation curves for Left: IC 2522 and Right: ESO 582-G12. The errors on 
the data points are typically <10% (Persic & Salucci 1995). The solid lines represent the rotation curve that best matches the observed 
shear at 10 kpc. For IC 2522 the dashed line represents a model rotation curve with a bulge that has an affective radius which is a factor 
of 2 larger than measured in the 1-D decomposition, in order to mimic the affect of a bar with a major-to-minor axis ratio of ~2. 



rate, favors a non-adiabatically contracted halo, with a 
low halo concentration and a massive, extended dark 
matter halo. The second, with a high shear rate, fa- 
vors an adiabatically contracted halo with a high NFW 
concentration, and relatively low virial mass. These two 
cases motivate the application of the shear-rate con- 
straint to a larger sample of galaxies. Since shear is 
related to the spiral arm pitch angle and spiral arms 
are clearly detectable in disk galaxies at z ~ 1 (e.g. 
Elmegreen, Elmgreen fc Hirst 2004), it is possible to 
constrain models of disk formation and dark matter halo 
structure back to a look-back time of ~ 7 Gyrs, and eval- 



uate how the mass distribution in disk galaxies changes 
with time. 

An important test of our modeling is how well the mod- 
eled rotation curves fit the data. Figure 9 shows the 
observed rotation curves from Persic & Salucci (1995) 
overlaid with our modeled rotation curves. In both cases 
the rotation curve is modeled well in the outer parts, 
where we are measuring the shear. In the inner regions, 
where the bulge is important, both rotation curve mod- 
els overestimate the observed rotation curve, although 
the model for ESO 582-G12 presents a much better case 
than that for IC 2522. For IC 2522, the solid line repre- 
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Fig. 10. — Ks band image of IC 2522, revealing a bar that was 
not eveident in the optical imaging of this galaxy. The image is 
aligned such that North is up and the field is 3.3x3.3 arcmin. 



sents a model rotation curve, with a bulgc-to-disk ratio 
of 0.164 (sec Tabic 4) and a bulge effective radius of 0.2 
kpc (as measured in the 1-D decomposition). However, 
infrared imaging of IC 2522 (see Figure 10) reveals that 
this galaxy has a bar with a major axis length ~2 kpc 
(i.e. the radius within which the model overestimates the 
observed rotational velocities) and a major-to-minor axis 
ratio of ~2. (It is common for near-IR imaging to reveal 
bars that may have been hidden due to dust extinction 
in an optical image, e.g. Seigar & James 1998a; Eskridge 
et al. 2000; Seigar 2002; Seigar et al. 2003). The spec- 
troscopic data taken by Mathewson et al. (1992) were 
observed with a slit aligned along the major axis of each 
galaxy, and in the case of IC 2522, the major axis of 
the bar is well aligned with the galaxy major axis. The 
stellar orbits within a bar are such that the dominant 
motion is parallel to its major-axis (e.g. Athanassoula 
1992), and so this will account (to some extent) for the 
low values of the measured rotation velocities (compared 
to the modeled values) within a few kpc. Furthermore, if 
we double the effective radius of the modeled bulge, the 
resulting rotation curve is the dashed line in Figure 9. 
This still fits the outer part of the galaxy rotation curve 
extremely well, and the central velocities are not as badly 
overestimated. Modifying the mass models in this way 
does not affect the trends displayed in Figures 6 and 8, it 
just affects how well the model reproduces the inner part 
of the observed rotation curve. This is probably because 
the main constraints we are using to model these galaxies 
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are at a 10 kpc radius, well outside the region where the 
bulge or bar is dominating the rotational velocities. As a 
result, the concentrations we determine from our model- 
ing do not depend on the mass distribution assumed for 
the bulge (i.e. the mass distribution within the central 
few kpc). This approach is, therefore, a powerful and 
robust method, because its constraints are relatively in- 
sensitive to the details of bars versus bulges in the central 
regions. 

6. CONCLUSIONS 

In this paper we have shown that near-infrared and op- 
tical spiral arm pitch angles are the same, on average, and 
as a result we have strengthened the correlation between 
spiral arm pitch angle and shear rate originally shown in 
Seigar et al. (2005) and expanded here using optical data. 
Using an infall model, we have shown that the use of rates 
of shear (which can now be derived from spiral arm pitch 
angles) allow us to put constraints on the total central 
mass concentration, the dark matter concentration and 
the initial NFW concentration. In some cases it may be 
possible to determine if the infall has to occur adiabat- 
ically or non-adiabatically, and this is demonstrated by 
IC 2522 which has to undergo non-adiabatic infall. 

This method can be used to determine the central con- 
centrations of galaxies and to constrain galaxy formation 
models in any galaxy that has detectable spiral struc- 
ture. In future papers we will apply our technique to a 
large sample of galaxies. Furthermore, since spiral struc- 
ture can clearly be seen in disk galaxies out to z ^ 1 
(Elmegreen et al. 2004), the evolution of central mass 
concentrations in disk galaxies can be estimated as a 
function of look-back time, and we plan to investigate 
this as well. 
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